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_( ■ Abstract. We present the first determination of the pairwise velocity dispersion of galaxies 

r- { ' at different luminosity with the final release of the Two-Degree Field Galaxy Redshift Survey 

(2dFGRS). Our result surprisingly shows that the random velocities of the faint galaxies are 
very high, around 700kms~^, reaching similar values as the brightest galaxies. At intermediate 
luminosities slightly brighter than the characteristic luminosity Af*, the velocities exhibit a 
well defined steep minimum near 400kms~^. The result challenges the current halo model of 
^ ■ galaxies of Yang et al. that was obtained by matching the clustering and luminosity function of 

1/^ ' 2dFGRS, and can be an important constraint in general on theories of galaxy formation, e.g., 

t~^ , the semi-analytical model. Combining the observed luminosity dependence of clustering, our 

^_^ • result implies that quite a fraction of faint galaxies are in massive halos of galaxy clusters as 

>0 I the brightest ones, but most of the A/* galaxies are in galactic halos. 

o 

:^ 



1. Introduction 



The clustering of galaxies in the Universe is characterized by their spatial positions, 

O ' and their peculiar velocities which lead to deviations of their motion from the pure 

Hubble flow. The big redshift surveys assembled in recent years by the diligent work of 

^ ' many astronomers give angular positions and redshifts for large numbers of galaxies. A 

rough 3D map can be obtained by placing the galaxies at distances along the line of 

sight derived via Hubble's law from their redshifts. The peculiar velocity, however, also 

^ , contributes to the redshift, and this leads to a misplacement of the galaxy away from its 

V^ ' true location. The local gravitational field is the cause of the peculiar motion, and thus 

the redshift distortion in the galaxy maps can give information on the underlying matter 

distribution. 

The amplitude of the distortions can be estimated from the pair distribution of galaxies. 
For pairs of galaxies at distances much larger than their separation, one can use a plane- 
parallel approximation, and write for the power spectrum in redshift space. 

Here /i is the cosine of the angle between the wave vector and the line of sight. The linear 
redshift distortion parameter (3 is related to the linear growth factor fiytM) — ^°*^ (^m 
is the matter density), and the bias b of the density fluctuation spectrum via (3 = Vt^^f /b. 
P(k) is the power spectrum in real space. The pairwise velocity dispersion (PVD) ay{k) 
describes the virial motion of galaxies in dense systems, i.e. the Finger-of-God effect. 
As shown by Jing & Borner (2001), the redshift power spectrum of dark matter and 
its biased tracers can be accurately described by Ea. (|l.l() . Applying this model to an 
observation of P^{k, /i) may therefore yield a determination of the three quantities P{k), 
/3, and ay (fc) , all of which are useful observables for testing galaxy formation models 
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Here we use the final release of the 2dFGRS ( http : //www . mso . anu . edu . au/2dFGRS) 
to measure P^{k, /i) and the PVD (Jy{k) of galaxies. Since we are interested in the PVD 
at small scale, we fix /3 at a reasonable value of 0.45 (Peacock et al. 2001). The 2dF 
catalog has already been analyzed statistically with respect to the PVD (Hawkins et 
al 2003), but here we use a novel method to estimate the PVD. Furthermore, we take 
advantage of the large number of galaxies in the 2dF catalog to bin the galaxies in 
different luminosity intervals, and to make the first study of the luminosity dependence 
of the PVD. This was not possible up to now, and we shall see that remarkable tests 
of the galaxy formation models become possible with the luminosity dependence of the 
PVD. The luminosity dependence of the clustering (the two-point correlation function) 
of galaxies in the 2dFGRS has been investigated (Norberg et al. 2002a), but not of the 
PVD. 



2. Statistical Analysis 

In order to study the luminosity dependence of the PVD, we divide the galaxies into 
11 subsamples according to their absolute luminosity. The subsamples are successively 
brightened by 0.5 magnitude from the faintest sample Mi, = —17.0 + 51og/i to Mt = 
—22.0 + 51og/i, with successive subsamples overlapping by 0.5 magnitude. Here h is the 
Hubble constant in units of lOOkms^^Mpc" . For computing the absolute magnitude, 
we have used the k-correction and luminosity evolution model of Norberg (2002b; k + e 
model);, i.e., the absolute magnitude is in the rest frame bj band at z = 0. We assume a 
cosmological model with the density parameter fio = 0.3 and the cosmological constant 
Ao = 0.7 throughout this paper. 

We measure the redshift two-point correlation functions £,z{s) following the method of 
Jing, Mo, & B6rner(1998). The random samples for the clustering analysis are generated 
in the same way as described in Jing & Borner (2004a). Each random sample for a 
northern or southern luminosity subsample contains 100,000 random points. We convert 
the Cz(s) to P^{k,ii) by the Fourier transformation. More details about our statistical 
method can be found in our journal paper (Jing & Borner 2004b). 

3. Results 

In Fig.l the basic measurement of the power spectrum in redshift space P(fc,/i) is 
shown. The four panels in this figure correspond to four different luminosity intervals 
from faint to bright galaxies. The values of k range from 0.16 /iMpc^^ at the top to 
2.5 /iMpc~ at the bottom with an interval A log^o k = 0.2. The south and north samples 
agree quite well with the full survey indicating that cosmic variance is not a major 
problem. The power spectrum for the larger k-values decreases quite strongly with /i, 
more than a factor of 10 between fi — and /i = 1. The solid lines are the best fits 
obtained by applying equation Hl.lfl to the data of the whole survey, showing that the 
model (|1.1() can accurately describe the observed redshift power spectrum of galaxies. 

The luminosity dependence of the PVD is shown most clearly in Fig. 2, where we 
have plotted ay at k = l/iMpc~^ for the 11 overlapping samples. The surprising result 
is a strong dependence on luminosity with the bright and the faint galaxies reaching 
high values of ~ 700kms~^ or more, and a well defined minimum of ~ 400kms~^ for 
the galaxies of magnitude Mi, — 51og]^Q/i = —20.5. The bright and the faint galaxies 
apparently have high random motions, as expected for objects in massive halos or in 
clusters. The M+ like galaxies are rather moderate in their PVD, and probably reside in 
galaxy size halos. 
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Figure 1. The redshift space power spectrum P{k,n) measured in 2dFGRS. Tlie symbols 
are for the whole survey, the dotted lines for the south subsample, and the dashed lines for 
the north subsample. The errors are plotted only for the whole survey, estimated with the 
bootstrap method. The smooth solid lines are the best fits of Eg. 11. Ill to data of the whole 
sample, (a) for -17.5 < Ah ~5logh < -18.5; (b) for -18.5 < Mt - 5 log ft < -19.5; (c) for 
-19.5 < Mb- 5logh < -20.5; (d) for -20.5 < Ah - 5 log /i < -21.5. The k values in each panel 
are from 0.16/iMpc~^ (top) to 2.5/iMpc~^ (bottom). 



The thick solid line represents the prediction based on the up-to-date halo model 
of Yang et al. (2003). For the halo model, we adopt the cosmological model that is 
a currently popular flat low-density model with the density parameter Jig = 0.3 and 
the cosmological constant Aq = 0.7 (LCDM). The shape parameter F — flgh and the 
amplitude as of the linear density power spectrum are 0.2 and 0.9 respectively. We use 
two sets of simulations, with boxsizes L = 100/i~^Mpc and L = 300/i~^Mpc, that were 
generated with our vectorized-parallel P^M code (Jing & Suto 2002) for this model. Both 
simulations use 512'^ particles, so the particle mass rup is 6.2 x IO^Mq and 1.7 x 10^'^Mq 
respectively in these two cases. We populate the halos with galaxies in a similar way to 
that of Yang et al. (2004), but we used our own code and adopted the model parameters 
of Model Ml in Yang et al. (2003). The halo model, which matched nicely the luminosity 
function and the luminosity dependence of clustering of galaxies in 2dFGRS (confirmed 
by our analysis), clearly does not match our observation of the PVD. The failure of the 
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Figure 2. The pairwise velocity dispersion measured at fc = 1 ftMpc" in the 2dFGRS (symbols 
for the whole sample, dotted line for the south, and dashed line for the north), compared with 
the prediction based on the halo model (thick solid line) . The error bars of the observed results 
are given by the mock samples. 

halo model indicates that the prescription of how to populate halos must be adapted 
better. 



4. Discussion 

The analysis of the velocity fields of the galaxies in the 2dFGRS has led to a surprising 
discovery: The random velocities of the faint galaxies are very high, around 700kms~^, 
reaching similar values as the bright galaxies. At intermediate luminosities the velocities 
exhibit a well defined steep minimum near 400 km s^^. 

It seems that the galaxies in different luminosity intervals appear as different pop- 
ulations in their own right, defined by objective statistics. A look at figure 2 shows 
convincingly that this is actually the case. For this figure we have sorted the galaxies in 
11 luminosity bins, each one magnitude wide, from magnitude —16.5 to —22 and plotted 
the value of a^ at a wave number of fc ~ 1 /iMpc^ . Such a finely resolved binning of 
galaxies in samples of different luminosities is possible for the 2dFGRS, because it is big 
enough to contain sufficiently many galaxies in each luminosity class. 
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The PVD is an indicator of the depth of the local gravitational potential. Therefore 
we find the interesting result that the bright and the faint galaxies move in the strongest 
gravitational field. They are in clusters, while the galaxies around magnitude —20, the 
M* galaxies in the Schechter luminosity function, populate the field. 

The bimodal nature of the correlation between the PVD and luminosity may be used 
as a stringent test of galaxy formation models. We have investigated the halo occupation 
model (Yang et al. 2003) which has been optimally fitted to reproduce the luminosity 
function, and the two-point correlation function of the 2dFGRS. If we adapt this model 
to the PVD value of the Af* galaxies, we see that it cannot give the high values found 
for the fainter galaxies. The PVD values of the model actually run opposite to the data 
and the model assigns smaller values to the fainter galaxies. This must mean that the 
assignment of galaxies to the dark matter halos must be done in a more intricate way as 
up to now. The number of faint galaxies in clusters must be increased substantially to 
at least recover the high PVD found for them. Also, the low value of 400kms~^ found 
for the galaxies with magnitude —20.5 must mean that these galaxies reside in dark 
matter halos of galactic size, The halo population model must be fitted with a much 
more complex scheme of assigning galaxies to halos, if the results shown in figure 2 are 
to be reproduced, 

Another way, widely used, to connect dark matter to galaxies is the semianalytic 
modeling, where the dark matter distributions obtained from N-body simulations are 
supplemented with some of the physical processes important in galaxy formation using 
semianalytic techniques. A test of the PVD vs luminosity for this type of models will be 
the aim of a subsequent paper. 

Acknovirledgements 

JYP would like to thank the Max-Planck Institute fiir Astrophysik for its warm hos- 
pitality, lAU for its partial travel support, and the SOC for adjusting the talk time. The 
work is supported in part by NKBRSF (G19990754), by NSFC (No.l0125314), and by 
the CAS-MPG exchange program. 



References 

Hawkins, E., et al. 2003, MNRAS, 346, 78 
Jing, Y. P. & Borner, G. 2001, ApJ, 547, 545 

& Borner, G. 2004a, ApJ (in press)/ |astro-pli/0311585l 

& Borner, G. 2004b, in preparation 

Mo, H. J., & Borner, G. 1998, ApJ, 494, 1 

& Suto, Y. 2002, ApJ, 574, 538 
Norberg, P. et al. 2002a, MNRAS, 332, 827 
Norberg, P. et al. 2002b, MNRAS, 336, 907 
Yang, X. et al. 2004, MNRAS (in press) /astro-ph/0303524 
Yang, X., Mo, H. J., & van den Bosch, F. C. 2003, MNRAS, 339, 1057 
Peacock, J. A., et al. 2001, Nature, 410, 169 



Jing, 


Y. 


P. 


Jing, 


Y. 


P. 


Jing, 


Y. 


P. 


Jing, 


Y. 


P. 



